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ABSTRACT 


An analytical investigation was made of a transient method 
of calibrating fine wire thermocouples. The system consisted 
of a 10 mil dia. standard thermocouple (Pt, Pt-13% Rh) and an 
0.8 mil dia. chromel-alumel thermocouple attached to a 20 mil dia. 
electrically heated platinum wire. The calibration procedure 
consisted of electrically heating the wire to approximately 2500°P 
within about a seven-second period in an environment approximating 
atmospheric conditions at 120,000 feet. Rapid periodic readout of 
the standard and fine wire thermocouple signals permitted a 
comparison of the two temperature indications. An analysis was 
performed which indicated that the temperature distortion at the 
heater wire produced by the thermocouple junctions appears to be 
of negligible magnitude. Consequently, the calibration technique 
appears to be basically sound , although several practical changes 
which appear desirable are presented and discussed. It appears 
that some additonal investigation is warranted to evaluate 
radiation effects and transient response characteristics. 



INTRODUCTION 


The most common method of measuring the temperature of a wall 
is to attach a small beaded thermocouple to the surface. In the 
present application it was desired to make such surfaces tempera- 
ture measurements with very fine wire thermocouples, on the order 
of 0.001" diameter. The actual application contemplates the 
measurement of the transient surface temperatures of re-entry 
spacecraft, such that the temperatures experienced range up to 
2500°F, while the environmental conditions of pressure and gas 
compositions are those of the atmosphere at re-entry at about 
120,000 feet. 

In order to insure the accuracy of thermal sensors utilized 
by the re-entry spacecraft vehicles it was desired to investigate 
the performance characteristics and accuracy of the fine wire 
thermocouple to be used. To accomplish this a test procedure 
was devised in which a 1-mil fine wire chromel-alumel thermo- 
couple, together with a 10-mil platinum-platinum-1 3% rhodium 
reference thermocouple, was attached to the center of a 20-mil 
platinum-13% rhodium electrically heated wire. The two fine wire 
thermocouple leads and one of the leads of the reference thermo- 
couple formed a single junction as shown in Figure 1. The test 
assembly was located within a vacuum chamber that permitted 
operation under vacuum conditions in various gaseous environments, 
as indicated in Figure 2. 

The calibration procedure consisted of rapidly heating the 
central wire to about 2500°F within about a seven-second period 
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and comparing the output of the two thermocouples. A transient 
testing procedure was considered essential since extended 
periods of operation at elevated temperature in an oxidizing 
atmosphere result in deterioration and damage of the thermo- 
couple materials. Because a transient technique of calibration 
was required it became imperative to perform an evaluation of 
the validity and accuracy of the results obtainable from such a 
technique. Specifically, the goals were to: 

1. Perform an analytical investigation of the heat 
transfer behavior and characteristics of the 
thermocouple heater wire junction. 

2. Estimate the errors to be expected. 

3. Present recommendations for improving the accuracy 
of the technique. 

The following sections will describe the efforts made toward the 
preceding goals. 



APPARATUS AND PROCEDURE 


The essential features of the transient calibration technique 
were presented in the previous section. Only additional supple- 
mentary information will be presented here. 

The thermocouple junction and associated apparatus are 
presented in Figures 1 and 2. The test procedure consisted simply 
of applying a pre-set current to the heater wire and obtaining 
simultaneous readings of the standard and test thermocouples at 
reading rates up to 40 measurements per second over a seven-second 
interval . 


ANALYSIS 

One of the first elements to consider in investigating the 
accuracy of any temperature measurement is the possible disturbance 
of the temperature field being measured produced by the measuring 
device itself. A second factor to consider here, since the 
calibration procedure is transient in character, is the response 
of the thermocouple to the rapid temperature variation experienced 
by the heater wire. 

Concerning the former factor, the possibility exists that the 
temperature field in the heater wire may be distorted appreciably 
in the region of the standard and fine-wire thermocouple junctions. 
Such a complex system of intersecting cylinders does not appear to 
lend itself readily to analytic solution. Instead, it appears 
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initially that the governing partial differential equations, 
together with their initial- and boundary-conditions, which 
could yield temperature -time-space data for the junction, is 
most amenable to solution by finite-difference approximations 
and numerical techniques. Initial efforts indicated that 
obtaining solutions by numerical techniques involved considerable 
computational difficulties (which are not insurmountable if time 
is not a factor) . These efforts are described in Appendix I. 

It was decided instead to obtain closed form analytical solutions 
of the temperature field for a simplified junction approximating 
the actual one. 

The simplified junction considered was that of a 0.010" dia. 
platinum wire attached to the 0.020" dia. heater wire. Since 
the cross-sectional area of the reference junction leads are over 
100 times larger than the fine wire thermocouple leads, it was 
believed that disturbance of the temperature field due to the 
reference lead would be much larger than that due to the test 
lead. Consequently, the presence of the test thermocouple leads 
was not considered initially for purposes of analytical tractability . 

The system considered, then, is illustrated in Figure 3. The 
governing equations for the two domain systems are: 
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With reference to equations (1) and (2) several observations 
should be made. Heat flow in the heater wire is assumed to be 
one-dimensional in the radial direction. This is a good assumption 
for wires of large L/d. Heat flow in the reference thermocouple 
is also assumed one-dimensional in the axial direction. Since 
the thermocouple has only one-fourth the cross-sectional area as 
the heater wire and behaves as a fin attached to the heater 
wire this assumption is believed justified. The implicit assump- 
tion is made that the thermocouple heater wire junction is flat, 
which is a reasonable approximation due to the relative sizes of 
the wires. 

The solutions to equations (1) and (2) are developed in 
Appendix II and are: 

Heater Wire Surface Temperature: 
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Reference Thermocouple Temperature: 


T-T* 




• X 


(5) 


T- T 

/X 

The solutions (3) and (5) were evaluated at 0.5 second 

intervals for a seven-second period and are listed in Figure 4. 

For a given instant of time the heater wire surface temperature 

served as one of the boundary conditions for the thermocouple 

(i.e., T (0,t)=T (R/t)). The temperature distribution in the 
1 2 

thermocouple was evaluated from (5) , which actually describes 
the steady-state distribution in a fin. Consequently, this 
procedure contains an implied assumption that the temperature 
throughout the fin reacts instantaneously to changes in the 
heater wire surface temperature. 

It is conceivable, however, that under actual test conditions 
thermocouple conduction, convection and radiation losses may drain 
sufficient energy from the heater wire so as to result in a lower 
base temperature than calculated from (4) and (5). Practically 
speaking, the temperature measured by the thermocouple would be 
lower than the undisturbed surface temperature of the heater wire. 
In order to estimate the possible difference between the 
theoretically determined temperature and the actual temperature 
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the following procedure was followed. As noted in Figure 4, the 
temperature experienced by the thermocouple node in contact with 
the heater wire after 3.5 seconds (or midway through the test) 
was 1317°F. Using this temperature and the steady-state tempera- 
ture gradient in the thermocouple at that time , the convective , 
conductive and radiation heat losses from the node were evaluated, 
in addition to the stored energy within the node. The effect of 
these heat losses upon the temperature distribution in the 
thermocouple lead was then determined for a seven-second test 
period. The modified temperature distribution in the test for a 
seven-second period are listed in Figure 5. Comparing the 
temperature distributions of Figures 4 and 5, it is seen that 
modification for conduction and radiation losses produce very 
little effect on the junction node temperature — a very important 
occurrence that will be discussed more thoroughly in a later 
section. 

Since the temperature experienced in these calibration tests 
range from 70 to 2500°F, accurate solutions to equations (1) and 
(2) require accurate thermophysical properties data and heat 
transfer coefficients over the entire temperature range. The 
thermophysical properties employed are presented in Appendix III. 
In general, the temperature dependence of all thermophysical 
properties was taken into account in obtaining transient tempera- 
ture histories from equations (1) and (2) . 
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The free convection heat transfer coefficient was initially 
evaluated from correlations for laminar and turbulent free 
convection from horizontal and vertical cylinders (Equations 
7-18, 7-21, 7-29, 7-30 from [l])for the heater wire and thermo- 
couple respectively. The coefficients so evaluated are not 
accurate, however, for the small wire sizes and high temperatures 
experienced here. A more accurate correlation is that of McAdams 
[2] , in which data are well correlated over a range of Grashof 
numbers from 10 - ^ to 10^ for fluids as different as air and 
glycerin, for cylinders as small as 0.003 cm in diameter. Further 
information regarding convective heat transfer coefficients is 
presented in Appendix III. 


DISCUSSION 

Temperature Distributions in Junction Region 

As indicated earlier, the analytically determined junction 

node temperature of the 0.020" reference wire is influenced very 

little by the average heat losses occurring at a temperature of 
o 

about 1300 F for the entire seven-second test period. This 
occurrence appears to imply that the temperature within the heater 
wire in the junetion region is not changed appreciably by the 
presence of the reference thermocouple. The results of Figure 5 
indicate that the conduction, convection and radiation losses 
from the reference junction node produce a temperature change of 
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only a small fraction of a degree over a seven-second period. 

Viewed in a different perspective, the energy generated within 
the heater wire within a 3.5 second period is of such a magnitude 
as to raise the temperature of the thermocouple approximately 
1245°F. In comparison the temperature change occurring as a 
result of the maximum heat losses within the reference junction 
node acting over twice that length of time produces only a very 
small temperature change. Consequently, the energy "drained'' 
from the heater wire by the reference thermocouple is of such 
comparatively small magnitude that the temperature distribution 
within the heater wire is not disturbed significantly by the 
presence of the reference thermocouple leads. If this is true, 
then the disturbance produced by the pressure of the fine-wire 
thermocouple leads is even less significant, since they offer a 
conduction path less than 1/100 times as large of that of the 
reference beads. Other possible sources of error, however, must 
be evaluated, one of which is the error associated with the precise 
nature of attachment to the heater surface. 


Method of Attachment 

The temperature gradient in the attached thermocouples, as 
exemplified by the temperatures tabulated in Figure 4, is quite 
steep close to the heater wire. This occurrence is related to 
another type of eeror that may exist in thermocouple measurements, 
namely that associated with the presence of a third material 
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forming the thermocouple bead or providing the means of attachment 
to a surface. Generally the possible existence of a temperature 
difference across a thermocouple bead or substrate material is 
neglected in theoretical models of thermocouple systems. That is, 
most analyses deal with the case where the thermocouple bead is 
very small. Such simplifications are not always justified, as 
indicated from several investigations. Suh and Tsai [3] , for 
example, found that in a typical measurement of a burning solid 
propellant the temperature difference across the thermocouple bead 
is of the order of 70 to 100°C, for a thermocouple of h mil dia- 
meter . 

Errors associated with thermocouple response characteristics 
also may be influenced by the presence of substrate or bead 
material. Such errors are discussed in [4], where tests were 
conducted on thermocouples surrounded with silver epoxy, plain 
epoxy, no epoxy, solder and glue- Not surprisingly, various 
response characteristics were observed. Perhaps of even greater 
importance relative to the present investigation is the 
following statement by Parker [4]: "In addition to the experimental 
work, several analytical heat transfer solutions were made with 
the use of a computer and compared with the experimental data. 

In general, there was poor agreement because it. was extremely 
difficult to determine very small thicknesses between the thermo- 
couple and substrate accurately." 


N 
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In the calibration technique investigated here, very high 
temperature levels and large temperature gradients near the 
junction region are experienced, in addition to rapid temporal 
variations. Consequently the preceding comments are quite 
relevant. In general, any uncertainty in the precise location 
of the thermocouple junction with respect to the heater wire 
surface, whether produced by the presence of substrate material 
or a welding process, may result in significant error. If a 
bead is formed at the junction, obviously the uncertainty 
associated with response characteristics may be cause for 
considerable error. 


Nature of Thermocouple Junction 

The most common type of thermocouple junction employed is 
the bead junction, wherein the temperature indicated is a function 
of the temperatures where the wires leave the bead. Practically, 
this translates into making the bead small and flattened, being 
in intimate contact with the surface whose temperature is being 
measured . 

In the present application the separated open junction is 
employed [5], That is, each wire is separately joined to the 
surface, which must be an electrical conductor. This type of - 
junction, which is actually two series junctions, delivers an 
output that is a weighted mean value of the two individual 
junction temperatures: 
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From equation (6) , the difference between the output and 
mean emf may be minimized by reducing the temperature difference, 
T'-T". This type of junction has been shown to be more accurate 
than the bead junction if this condition is adequately fulfilled 
[6,7] . In the present case, due to symmetry, the reference leads 
should be at essentially the same temperature , if they are truly 
surface welded and no substrate material (or differences in 
substrate material) is present at the junctions. 

It appears, however, that placement of the fine wire 
junctions is not established as precisely. In particular, it 
appears quite possible that one junction may be totally in contact 
with the heater wire surface immediately adjacent to the reference 
junction, while the second may be in partial contact with the 
reference thermocouple wire. Thus, since T' jt T" in both cases, 
the indicated temperatures of the two thermocouples would always 
differ, even if all other errors we re either eliminated or 
negligible. A more judicious placement of the fine wire thermo- 
couple leads appears warranted. 
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Response Characteristics of Thermocouples 

Both analytical and experimental techniques have been 
developed in the literature to describe the transient response 
characteristics of thermocouples, [8,9]. In the present 
application the accuracy of the calibration technique appears 
strongly dependent upon the comparative response characteristics and 
may be of greater relevance than errors due to conduction and 
radiation losses per se. In the present calibration procedure 
the possibility exists that the test thermocouple responds more 
quickly to the heater wire surface temperature than the relatively 
more massive reference thermocouple. Thus, though a difference 
between the indicated temperatures of the two thermocouples may 
be observed, it is conceivable that the temperature of the fine 
wire thermocouple may be more accurate than that of the reference 
thermocouple which was initially intended to provide the "true" 
surface temperature. 

Generally, Laplace transform techniques are used to analyze 
the response characteristics of thermocouples. The indicated 
temperature and actual temperature are related by the transfer 
function as 


G (s) 

m 



(7) 
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T. 


M 




= thermocouple transfer function 

= Laplace transform of the measured temperature 

= Laplace transform of the actual or calculated 
temperature 
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For a thermocouple welded to the external surface of an electrically 
heated test-section G T (s) is of the form {9] 


G(s) * —S— 
T s + c 


( 8 ) 


where a and c depend upon the thermocouple installation and may 
differ for each thermocouple. It is possible to experimentally 
determine these values for particular thermocouples by imposing 
a known temperature change (ramp change, for example) upon the 
test section; from the known values of the indicated and actual 
(calculated) temperatures, together with their Laplace transforms, 
it is possible to calculate the values of a and c in (2) , thus 
determining G T (s). Once the transfer function has been established, 
it becomes possible to calculate the true or actual temperatures 
from the indicated or measured thermocouple temperatures from the 
relation 


T 



A 

TsT 


(9) 


The present investigation has bean concerned primarily with 
the transient temperature distributions in the heater wire and 
thermocouple junction region, and the related inaccuracies due to 
conduction and radiation losses. It is this writer's judgment that 
additional time should be devoted to investigation of the 
response characteristics of the thermocouple assembly before the 
accuracy of the calibration technique may be precisely and 
quantitatively evaluated. 
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Radiation Effects 


If radiation heat losses from the reference thermocouple and 
heater wire are taken into consideration, equations (1) and (2) 
are modified as follows 


Heater Wire: 
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It is seen that equations (10) and (11) are both non-linear, 
(10) by virtue of the radiation boundary condition, and (11) 
by virtue of the non-linear radiation term in the governing 
equation itself. Solution of this system of non-linear partial 
differential equations is a considerable level of difficulty 
higher than that required in the solution of (1) and (2) , and 
must be done numerically. 
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Even if (10) and (11) were to be solved/ however, a factor 
which is not taken into account by these equations is the 
radiation exchange between the heater wire and thermocouple lead(s). 
Inclusion of this closer approximation to reality involves an 
additional considerable increase in level of difficulty in solving 
the governing system of equations. In view of these inadequacies 
of analytical techniques it appears simpler, and advisable in this 
writer's opinion to evaluate radiation effects by suitably devised 
experiments. 


Deterioration of Thermocouple Performance at High Temperatures 

Initial attempts to calibrate fine wire thermocouples at 
temperatures above 1220°F proved difficult since a period of 
about three minutes was required to obtain each set of tempera- 
ture measurements. The oxidation occurring at these high 
temperatures produced deterioration of the thermocouples, thereby 
producing the need for the rapid transient calibration procedures 
being investigated here. Although the time at which the thermo- 
couples are at high temperatures is short, the possible influence 
of surface chemical reactions upon thermocouple accuracy should be 
considered. 

Reference [10] discusses at some length catalytic activity 
in platinum group temperature sensors in connection with possible 
measurement errors. As in several other investigations the 
authors proposed the use of non-catalytic coatings on thermocouples 




17 


for most reacting gas temperature measurements. It appears 
advisable in the present transient calibration procedure to 
also employ such coatings. In addition to preventing or 
minimizing chemical decomposition of the thermocouple wires , 
the radiation heat interactions between the heater wire, 
thermocouple wires and environment, which are quite difficult 
to accurately determine analytically, would be greatly reduced 
due to the lowered surface temperature of the coating. The 
convective heat losses from the thermocouple leads would also be 
somewhat reduced for the same reason. Addition of coating 
material, however, could affect the relative response character- 
istics of the two thermocouples. The ultimate effect upon the 
accuracy of the calibration technique produced by the addition 
of coatings appears to be primarily dependent upon the two 
pnenonmena — reduced radiation heat transfer losses and transient 
response characteristics (although conduction losses and 
convective heat transfer with the environment must also be 
considered). Again, the need for investigation of transient 
response characteristics is indicated. Perhaps it would be 
possible to sufficiently minimize radiation losses from both 
thermocouples with a minimum coating thickness; further addition of 
coating material to the fine wire thermocouple might then be 
made in such an amount as to match the response characteristics 
of both thermocouples, resulting in an accurate calibration 
technique over the entire range of temperatures. 
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As to the particular type(s) of coating that should be 
employed, several substances have been employed as high 
temperature coatings in the recent literature, among which are 
quartz, silica, tetraethyl lead, molten borax, magnesia, 
oxidized chromium (C^Cl^) , NBS Ceramic Coating A-418, and 
others til] . The particular coating employed must of course be 
capable of withstanding the temperature range to be experienced, 
should have a low thermal conductivity to (a) minimize conduction 
losses, (b) to decrease the surface temperature and (c) to 
minimize aging effects caused by adsorption and desorption 
processes and exposure to a variety of atmospheres, and should 
of course be readily depositable upon the thermocouples. A 
more extensive discussion of thermocouple coatings and a compila- 
tion of pertinent literature may be found in [10] . 


CONCLUSIONS 

Although a detailed mapping of the temperature-time-space 
field within the heater wire thermocouple junction region has 
not been presented herein the analysis performed indicated that 
the temperature depression caused by the reference thermocouple 
leads being attached to the heater wire is negligibly small 
throughout the test. Consequently, the test thermocouple leads, 
having a cross-sectional area 1/100 times as large as the 
reference thermocouple leads, will have essentially no effect upon 
lowering the heater wire temperature. 
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Errors associated with the presence of substrate or bead 
material may be quite significant, due to the existence of a 
large temperature difference across the material or due to 
changes in response characteristics. Uncertainties in the 
placement of the fine wire junctions may also result. in error, 
indicating a more judicious placement of these junctions appears 
warranted . 

Errors due to a difference in response characteristics of 
the reference and test thermocouples may be more significant than 
that due to disturbance of the heater wire temperature produced 
by the presence (attachment) of the thermocouple wires to it. 

It is possible to experimentally determine the transfer function 
of each thermocouple. Additional experimental and theoretical 
investigation of the comparative transient response characteristics 
of the thermocouples also appears warranted. 

The addition of high temperature chemically inert coating 
materials to the thermocouples appears desirable principally in 
order to minimize radiation heat transfer interactions at high 
temperature levels and to minimize aging effects produced by 
exposure of bare thermocouple wires to various atmospheres. 


RECOMMENDATION S 

For reasons presented herein it is recommended that the fine 


wire junctions not be placed in the vicinity of the reference 
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junction (s). Since the reference thermocouple apparently disturbs 
the heater wire temperature very little, the fine wire leads need 
not be far removed. A distance of one inch should be more than 
ample. Both leads in each system should be welded to the surface 
diametrically opposite each other. Since the heater wire is seven 
inches in length it appears that conduction losses to the terminals 
should not influence the wire temperature near its center, where 
the two thermocouple systems will be located (more precisely, each 
thermocouple system will be located 2% inches from the voltage 
terminals) . The temperature depression in the heater wire produced 
by conduction losses to the voltage terminals may be easily 
calculated using expressions found in (12,13,14). 

Both sets of thermocouple wires should be coated with 
chemically inert high temperature coating materials to minimize 
radiation heat losses and to minimize deterioration and aging 
effects of the thermocouple wires. Addition of such materials, 
however, influences the transient response of thermocouple systems. 
Whether or not the thermocouples are coated the comparative 
transient response characteristics of the thermocouple systems 
should and must be analyzed before a quantitative evaluation of 
the transient calibration technique may be realized. In fact, 
the choice of a particular coating material may be a strong 
function of its influence upon the transient response character- 
istics of the thermocouple systems. 

The preceding recommendations have been analytically based 
or oriented. However, it appears essential to this writer that 



21 


an experimental program planned co-currently with the previously 
recommended activities be conducted. A wealth of information 
should be obtainable from a modest experimental program. For 
example, the transient response characteristics of the thermocouple 
systems with different coatings and coating thicknesses, indicated 
by analyses, could be experimentally evaluated. 

Likewise, the radiation heat transfer characteristics, which 
are rather difficult to establish quantitatively and accurately 
(due to mathematical complexity and lack of radiation properties 
data), could also be investigated experimentally. In addition, the 
comparative effects of radiation and convection losses could be 
established by performing tests under varying degrees of vacuum. 
Also, aging and deterioration effects may be possible by judicious 
testing in various gaseous enviornments at varying degrees of 
vacuum . 





Figure 1. Schematic of heater wire-thermocouple 
junction and instrumentation. 




Figure 2. Schematic of experimental apparatus 
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Figure 4 


Temperature distribution in thermocouple wire for 
seven-second heating period. 
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Figure 5 


Modified temperature distribution in 
thermocouple wire due to conduction and 
radiation losses. 
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NOMENCLATURE 

cross-sectional area of thermocouple 

slope of temperature-emf relation for two different 
materials, i.e. e = a Q + b T 

B iot number, h L/h 

thermocouple output if both junctions were at the 
mean temperature, (T' + T'')/2 

thermocouple output if junction of wire 1 and 
surface is at T' and that of wire 2 is at T" 

convective heat transfer coefficient 

thermal conductivity 

m = (hP/kA) 35 

perimeter of thermocouple cross-section 
radial coordinate 
radius of thermocouple 
Laplace transform variable 
temperature 

analytically evaluated surface temperature of 
thermocouple 

environment temperature 

temperatures of thermocouple wire-surface junctions 
in an open junction thermocouple 

internal heat generation rate per unit volume 

axial coordinate 

thermal diffusivity 

eigenvalue 

temperature difference, T-T^, 
time 
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Subscripts 

1 reference thermocouple lead wire 

2 heater wire 



APPENDIX A. 


Numerical Solution of Temperature Field of Junction Region 

Most of the material presented in this section is of a 
negative nature. That is, the approaches described here did not 
lead to the desired temperature -time-space information for the 
junction regions and are presented here for the sake of complete- 
ness. 

Initially it appeared that the most direct method of deter- 
mining the temperature history within the geometrically complex 
junction region, with non-linear (radiative) boundary conditions, 
and with temperature dependent thermophysical properties, was by 
means of numerical techniques. This contention remains valid. 
However, sufficient practical and even mundane difficulties were 
encountered in attempting computer solution of the governing 
equation as to make such solutions unattainable within the time 

i 

limitations set by the terms of this study. Some of the problems 
encountered will be described here. 

A cursory inspection of the geometry of the junction seems to 
indicate description of a nodal system in terms of cylindrical 
coordinates. It quickly becomes apparent that such a system is 
not workable because of extremely complex "communication" of 
nodes at the intersection of the thermocouples and heater wire. 
Consequently, a cartesian coordinate system was chosen in which 
the circular cross-section of the wires was approximated by a 
number of small square areas. 
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A system v/as chosen that compromised the numerical complex- 
ities of a very finely graded network and the approximate 
temperature distribution obtainable by using a very coarse 
network. A total of 155 equaions was developed describing 
the nodal heat balances and the boundary conditions. The 
equations were written in implicit form. This is necessary 
since the stability criteria derived considering the small wire 

sizes involved dictated solution of the 155 equations in time 

-5 

increments of roughly 10 seconds. For the entire test 
period of seven seconds, then, a prohibitive amount of computa- 
tions was required. The implicit method on the other hand is 
not restricted by the size of the time increment, but does 
involve the simultaneous solution of the 155 equations. 

A computer program using the subroutine SIMEQ from the 
Langley ACD library was developed to solve the system of 
equations. Due to inherent inaccuracies associated with the 
method of solution used in the subroutine it was not possible 
to obtain accurate nodal point temperature histories. Conse- 
quently a Gauss-Seidel iteration scheme, which has proved to 
be an accurate method of solution for even larger systems of 
equations, was employed. However, unexpectedly, problems of 
stability, associated with a parameter appearing in the 
governing equation, arose. At this point it was believed that 
time did not permit further investigation of the numerical 
approach. Thus it was decided to apply analytical techniques 
to an idealized model of the system to obtain results that would 
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provide some means for evaluating the accuracy of the transient 
calibration technique. 
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APPENDIX B. 


Analytical Solution of Heater Wire and Reference 
Thermocouple Temperatures 


The question which the following analysis attempts to answer 
is "Does the attachment of the reference thermocouple lead to the 
electrically heated wire provide a large enough path for heat 
leakage as to significantly depress the heater wire temperature 
over what the surface temperature would be without the thermocouple 
lead attached?" The simplified system, in which the presence and 
the effects of the fine wire thermocouple leads are neglected, is 
shown in Figure 3. The governing equations and boundary conditions 
are repeated here for convenience. It was assumed that the heater 
wire was infinite in both directions, and thus the temperature 
gradient in the axial direction was assumed zero. Also one- 
dimensional heat flow, in the axial direction, was assumed to exist 
in the smaller diameter reference wire. The governing equations 
are : 


Heater Wire: 



a 2 3T 



39 2 . u 1 * 1 
Tr k 2 


(B-l) 


I.C. s e 2 (r ,0) = 0 

B.C. 02(O*t) is finite or 

30 2 ( 0 , t ) * 0 

3r 


(B-la) 

(B-lb) 


(r,t) = h 2 e 2 


(B-lc) 



30 - 


Reference Thermocouple Lead: 


1 39 1 

3 2 0 

m 2 0 X 

(B-2) 

0^ Ft” 

3x 2 

X.C. : 

e 1 (x,o) 

= 0 

(B-2a) 


9 (0,T) 

= 0 (R, t) 

(B-2b) 


0 (°°,T) 

= 0 

(B-2c) 


To solve equation (B-l) and its boundary conditions the 
Laplace transform technique is conveniently employed. Taking the 
Laplace transform of (B-l) , 

Ee . aJ. 9 . ) = ail,, + k as + Hl l! 

a a 3 rz r 3 r k p ( B _ 3) 

where L (0 ) = O’ 

L(||) = p 0 - 6 (0) 

0(0) = 0 from (B-la) 

2 

Letting P/a = q , and rearranging (B-3) , 


r 


2 


d 2 0 d0 

dr 2 ” + r dr 


q 2 r 2 F = 


-u* * ' 
k p 


The solution to (B-4) may be written as 
The complementary solution 


F * 6 + 0 . 

c p 

is given by 


0 c = C i I Q (qr) + C 2 K Q (qr) 


(B-4) 


(B-5) 
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From (B-lb) , 0 is finite at r = 0, C 0 = 0. Thus 

C £ 


0 = C. I (qr) 

c 10 

2 

To find the particular solution, let 6 =» Ar + Br + C, 

P 

Substituting into (B-4) , one finds 


A = 0, B = 0, C = 


u 


1 I t 


u' " 


kq 2 p kp 2 


Combining the particular and complementary solutions. 


9 = C x X 0 (qr) + 


From (B-lc) 


-K|i= h a 


90 - 

or -k ^ = h 0 


@ r = R 


Applying the transformed boundary condition to (B-8), one 


c. 


Jk. f < 

1 


/;/ 


Substituting (B-10) into (B-9) and rearranging. 




-A. 




u'" 


) 


A 




(B-6) 

(B-7) 

(B-8) 

(B-9) 

obtains 

(B-10) 

(B-ll) 



One may now determine 9 by taking the inverse Laplace transform 


of (B-ll) , 


Q, 


Jk 



(B-12 ) 


The integral of (B-12) is seen to have a second order pole 
at z = 0 and a simple pole 


h X 0 {(|) H R> + Mf-)* 2 I 1 {(|)^R> * 0 


(B-13) 


Letting z = -aX n 2 / » one obtains 


J (X R) - X R J (X R) = 0 (B-14 ) 

k On n In 

which is a multiple value equation solved by many eigenvalues 

of X . 
n 

Now one may find the residues of the integral of (B-12) . 

Thus 

TZ 

h I (qr) e 

[ 2 ] (B-15) 

h l Q (qr) + hq I^(qr) 


Res f(z) = lim t 
at z = 0 z-*0 


2 

where z = P = q a 
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~f( 3 J __ j. L 

-O ~ 



and differentiating# 

H y) = j >A 7 a (i//x) + -kfJ, (%■/,<)( 

6 f° » , tL\ S 

l Aa I t (f4)e**X A J - 

Az* (^^)e a ^[izAi l CfR) + fA/? T 0 (% /<)j 

[A To (f-X) +- £ ^ T . (f/X] % 


Now as g 


Thus 


I (qr) * I (qR) = 1 
0 0 

(qr ) = (qR) = 0 


Res 

f ( 2 ) = lim 

q=0 

q 0 

Res 

f ( 2 ) = — — 

q=o 

2a : 


1 

2aq 


r- 


Now, to find the residue at z = - aX^ , for simplicity 

_ a x ^ ^ 

P = h I^(qr)e n 


Then 


Q = h I n (qR) + k q I (qR) 
0 1 


_ .. \ > 


/) < 

J6Uyr\ 




fl (- a. 


C ) 


£- 




* v /\ 




cv-rj/ 


■*_«)L tf> Tii 7 

. // y J 

l iMJdil e ~' N ~ T 

** a: [a u, ( f *)+*i 

r 


j A 1 ((- ~.*)e 


A Kl r 


«- * a .„ 3 Uzn^e)* d To ( >w * j _7 


(B-16) 


(B-17 ) 


(B-18) 


(B-19) 


let 


(B-20) 


(B-21) 
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But 


hT? 

X R J (X R) = ~ J (X R) 


n 1 n 


0 n 


and letting 


Bi = 


hR 

k' 




t; < h/yy 'T~ 


3 -"* 1 KLe; * a /? »J x 


(E 


And finally. 




s 


a"'«. 


^ e. 

— <• T~ 


V ^ 

- a \ A_. <• 




a 


J. 


/>1 




In terms of temperature, T, at any position within the heater 



<jO 


, /"r- 

T « 7 ~ +. ^ L r * 
/Z A ^ ^ — 




"'/? _ m'"6 ,'\ 

x A 


A 
e 


^7* 
✓vx C 


/7s ^ / 


<B 


In dimensionless form, 

T - T 

« _ zr , _L_ 

/-< ** a e: 

~ A 


-o 


v i 




•>< a r 

/h 


Z. 


sr< -/ 




(B 


- 22 ) 


-23) 


wire. 


-24) 


-25) 


-26) 
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For the solutions (B-23-26) , the values X R are zeroes of 

n 

J (X R) = 0 and are 
o n 

X x R = 2.4048 
X 2 R * 5.5201 
X 3 R - 8.6537 

X. R =11.7915 

4 

X R =14.9309 

5 


Treating the reference thermocouple lead as an infinitely 
long fin of circular cross-section/ the governing equation and 


its boundary conditions are: 

Hf.' _ yrrf&. ^ JL If. 


(B-27) 


where 


6 = t - T„ 


I.C. : 


6 1 (x / 0) = 0 

8, (0,t) = 6 (R,t) 
1 2 

8^(<»/t) = 0 


Assuming a product solution, 6^ = X(x)^(t), one obtains the 
ordinary differential equations 

X 1 ~ Y 2 at = 0 (B-28) 

where Y 2 = (8 2 -m 2 ), 6 2 is a separation constant and 

m 2 = hP/kA, 

and X" - m 2 X = 0 (B-29) 

whose solutions are, respectively. 



a y 2 t 

T <T> = e 
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Since 
Thus , 


mx -rax 
X (x) = Ae + Be 


6 , approaches zero for large values of X, A 
1 

-mx a y 2 t a Y 2 T-mx 

0 = B e e = B e 

1 


(B-30) 


(B-31) 


0. 


(B-32) 


Now at x * 0, 0=0, so that 

o 

a y 2 t 

0 = Be (B-33 ) 

0 


But 0 is also given by (B-23) evaluated at r = R, i.e. 0 (R) , 
o 2 

or a Y 2 t 

0 (0,x) * Be =0 (R,t) 

e 2 


Substituting into (B-32) , 

-mx 

0 (x,t) = 0 (R/r) e 
1 2 

where m = y 

/ kA 

Thus 

p -mx 

T(x,t) = T + T (R,t)-T |e (B-34) 

00 t 2 00 i 


As indicated in the body of the report, this solution contains 
the tacit assumption that the temperature distribution within the 
reference lead reacts instantaneously to changes in surface 
temperature of the heater wire. This is obviously not possible 



physically. However, for purposes of estimating heat losses 
caused by the presence of the thermocouple, this idealization 
is not unreasonable. 
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APPENDIX C. 

PROPERTIES 

For solids experiencing large temperature variations the 
heat flow and temperature distributions within those solids can 
be significantly influenced by a temperature dependent thermal 
conductivity (to a lesser degree the same holds for density 
and specific heat). Likewise, for a solid subjected to radiation 
and natural convection boundary conditions, both the radiative 
characteristics of the surface and the free convection heat 
transfer coefficient will in general be temperature dependent. 
Quite often, thermophysical properties data, radiation surface 
characteristics and accurate convection correlations are not 
always available. This will become evident in the following 
sections. 

Thermophysical Properties 

Thermophysical properties for chromel-P and alumel were 
obtained from J16) for the temperature range from about 200 to 
1000°F. Thermal conductivity values are shown graphically in 
figure C-l. For temperatures beyond 1000°F extrapolated values 
were employed by using the straight line equations. indicated 
below and also shown in figure C-l. 

Alumel: k =9.83 (10 ^)T + 14.58 Btu/hr-ft-°F 

3 (T in °F) 

p = 537 lb m /ft J 

c = 0.125 Btu/lb -°F 
p m 
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-3 ° 

Chromel-P: k =5.96 (10 )T + 9.76 Btu/hr-ft- F 

p = 545 lb /ft 3 

m 

c =0.107 Btu/lb -°F 
p m 


For platinum and platinum-rhodium alloys property values may 
be found in {17] . For the system of units used herein the 
following expressions describe the property vaaiations, as 


recommended in [17] : 

Platinum: k = 41.137-0.1387 (10" 2 )T 

p = 1340 lb m /ft 3 

c_ = 0.371 (10" 5 )T + 0.0311 

lr 

ct = 0.864-0.2913 (10” 4 )T 


Btu/hr-ft-°F 
(T in °F) 


Btu/hr=f t-°F 
(T in °F) 

f t 2 /hr 
(T in °F ) 


Platinum-10% Rhodium: 

k = 41.3 [1 + 1.1K10" 4 ) (T-64 .4)3 Btu/hr-ft-°F 

o (T in °F) 

p = 1266 lb/ft 

c =0.0383 Btu/lb -°F 
p m 

a = 0.853 tl + 1. 11 (10 -4 ) (T-64. 4)3 ft 2 /hr 

(T in °F) 

Since properties were not available for Pt-13% Rhodium, the 
foregoing properties for Pt-10% rhodium were used. 


Emissivit^ 


Total and spectral emissivities of platinum may be found in 


[17]. Emissivity, as understood therein, referred to the ratio 
of energy radiated hemi spherically per unit area per unit time by 
a surface to that radiated by a black-body at the same temperature. 
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Expressions describing the total emissivity of platinum are: 

-4 

e = 0.828 (10 ) T 240 <T< 1690 

e = 0.0837 + 0.333 (10” 4 )T 1690 <T< 2890 

(T in °F) 

No similar data was available for platinum-rhodium alloys. In 
the absence of such data the above expressions could likely be 
used without appreciable error. 

No data describing radiative surface characteristics of 
chromel-alumel materials appears available in the literature. 

Resistivity 

Should further computer solutions to the governing heat 

transfer equations be required in which variations in electrical 

resistivity with temperature of the heater wire be taken into 

account in order to simulate the actual calibration test, such 

values may be found in tabulated or equation form in [171. 

For temperatures up to 1500°C, an equation relating the 

resistivity of platinum to temperature is: 

-3 -7 

R. * R [1 + 3.9788(10 )T - 5.88(10 )T] ohms/cmf 

t ° (T in °C) 


Free Convection Heat Transfer Coefficient 

Generalized correlations for free convection heat transfer 
coefficients involve the property groupings of the Prandtl and 
Grashof numbers of the surrounding fluid. These properties are 
most often evaluated at the film temperature, i.e. the average of 
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the heated surface temperature and the temperature of the ambient 
fluid. Consequently, large variations of surface temperature 
affect the film properties upon which the heat transfer coefficient 
is based. In the present case, since the magnitude of the heat 
transfer coefficient obviously influences the temperatures experi- 
enced by the heater-wire-thermocouple assembly, the properties 
of the surrounding medium must be evaluated throughout the transient 
heating period. 

For air the properties c , k, p and the property grouping 

P 

g0p 2 /y 2 are presented in Figures C-2 to C-5 . Also shown in each 
figure are straight line approximations to various segments to the 
curves for use in computerized calculations. 

Another method for calculating properties as temperature 
levels change is to employ the subroutine shown in Figure C-6. 

This subroutine, when called upon to calculate a property at a 
particular temperature searches the computer memory to find the 
two property values at the temperatures immediately bracketing 
the temperature in question. A linear interpolation between the 
two temperatures is then made to compute the unknown property 
value. Since many temperature-property values may be initially 
read into the computer memory this straight line approximation 
is quite accurate. With this subroutine it is thus not necessary 
to use equations of the type shown in Figures 2-5. 
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Another alternative for calculating the heat transfer 
coefficient is d function of the ckangihg tempdiratute is to 
use the applicable correlation (p. 17(5, [2] in this instance) 
to calculate 'h' for the different wire sizes at several tempera- 
tures. Plotting these values, as in Figure C-7, one may either 
employ the equations of the curves drawn through the calculated 
points or make useeof the subroutine previously described. 



THERMAL CONDUCTIVITY k BTU/HR-FT- 
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Figure C-2. Specific heat of air at atmospheric 

pressure as a function of temperature* 
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Figure C~6. Subroutine VALUE 
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function of temperature. 


